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Abstract: If SUSY is discovered at the LHC, the task will immediately turn to determin- 
ing the model of SUSY breaking. Here, we employ a Mixed Modulus-Anomaly Mediated 
SUSY Breaking (MMAMSB) model with very similar LHC phenomenology to the more 
conventionally studied Constrained Minimal SUSY Model (CMSSM) and minimal Anomaly 
Mediated SUSY Breaking (mAMSB) models. We then study whether the models can be 
distinguished and measured. If we only fit to the various mass edges and mass end-points 
from cascade decay chains that are normally studied, a unique determination and mea- 
surement of the model is problematic without substantial amounts of LHC data. However, 
if event rate information is included, we can quickly distinguish and measure the correct 
SUSY model and exclude alternatives. 
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1 Introduction 

The experiments at the Large Hadron Colhder (LHC) are currently running and have 
begun to explore physics at the Terascale. In particular, the LHC has already probed large 
regions of the parameter space of supersymmetry (SUSY), one of the most well motivated 
and thoroughly studied approaches to solving the hierarchy problem [1-13]. 

Generically, the dominant signature for R-parity-conserving SUSY at the LHC is 
events with high transverse momentum (pt) jets (and maybe leptons) accompanied by 
large amounts of missing transverse energy (EJp^^^), originating from the pair production of 
squarks and gluinos which then cascade decay eventually to the lightest neutralino (LSP). 
Indeed, the strongest limits to date on the supersymmetric parameter space come from 
searches for events with multiple hard jets plus missing energy and zero or one lepton 



If new physics consistent with SUSY is discovered at the LHC, these cascade decay 
events may also be instrumental in determining the Lagrangian parameters of the underly- 
ing model. Various kinematic edges can be constructed from the jets, leptons, and missing 
energy. The measurement of these edges then constrains the SUSY mass spectrum [14- 
17]. Combined with additional measurements of sparticle decay branching ratios, these 
measurements can be used as input into a global fit for the parameters of a given SUSY 
model. 



[5, 7, 12, 13]. 
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Several groups have performed such fits for various supersymmetric parameter spaces 
[18-23]. In a sufficiently low-dimensional parameter space such as the constrained min- 
imal supersymmetric standard model (CMSSM) and with sufficient luminosity, the LHC 
measurements are powerful enough to enable a precise determination of the best-fit SUSY 
parameters. With low luminosity and/or in a parameter space with a greater number of 
degrees of freedom, however, the endpoint and branching ratio measurements alone are not 
always sufficient to ensure a stable fit and a precise parameter determination. 

In Ref. [18], it was demonstrated that the inclusion of a cross section measurement — in 
particular the rate of high-missing-p^ events — can drastically improve the precision of a fit 
to a SUSY parameter space with non-universal gaugino masses. The authors of Ref. [18] 
use a standard leading-order Monte Carlo simulation and a fast detector simulation to 
compute the rate at each point in parameter space. While this method is straightforward 
and general, there are some technical challenges to this approach that have so far made it 
impractical to include rates in most such fits. 

The first issue is the computational cost of a Monte Carlo simulation. In Ref. [18], this 
is addressed by only generating a small number (1000) of events per parameter space point, 
and by implementing a parallelized version of the Monte Carlo generator Herwig [24] on 
multiple processors. Even so, the computation of rates is extremely time-consuming and it 
is only feasible to scan a relatively small number of parameter space points when doing the 
fit. Because of the limitation on the number of events that can be generated per point in 
parameter space, there are also large statistical fluctuations in the computed rate. These 
fluctuations can spoil the stability of the fit. Finally, while the state-of-the-art calculations 
of SUSY production cross sections at the LHC include many higher-order corrections, the 
Monte Carlo generators only use leading-order cross sections. 

Recently, a technique to include rate information in fits and overcome the issues listed 
above was presented in Ref. [25], and was implemented in the fitting code Fittino [26]. In 
this implementation, the cross section for sparticle production is first interpolated from a 
grid in the space of sparticle masses that was pre-computed using the program Prospino 
[27]. The remaining factor in the event rate, the acceptance — the fraction of sparticle 
production events whose decay products pass the kinematical cuts — is computed using a 
novel semi-analytical method. Because computing the rate in this way is extremely fast 
and sufficiently accurate, it allows event rates to be used on the same footing as the usual 
edge measurements in fits to LHC data. We use this technique as implemented in Fittino 
in the fits performed here and find that the rates are often a crucial ingredient in obtaining 
convergent and precise measurements of SUSY parameters. 

Of course, we do not a priori know the mechanism of supersymmetry breaking, and thus 
which SUSY parameter space is the appropriate one in which to run such a fit. Most fits 
have concentrated on the CMSSM parameter space, because it involves a small number of 
parameters and because it has been extremely well studied. One could of course ask, would 
LHC measurements provide sufficient information to select one model of SUSY breaking 
over another? In other words, would a fit performed in the "wrong" parameter space yield a 
significantly worse goodness-of-fit than one performed in the "right" parameter space? This 
question was recently investigated [23] for a number of popular SUSY breaking models: the 
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CMSSM, which is a version of gravity-mediated SUSY breaking; minimal gauge-mediated 
SUSY breaking (mGMSB); minimal anomaly mediation (mAMSB) and large- volume string 
compactification models. Assuming a signal with 1 fb~^ of LHC data at -y/s = 14 TeV 
it was found that the CMSSM scenario can be distinguished, using the LHC endpoint 
measurements only, from all other scenarios considered except for the mGMSB, which can 
fit the data considered equally well. 

One SUSY breaking scenario that until now has not been studied in the context of 
a global fit is mirage mediation, also known as mixed modulus-anomaly mediated SUSY 
breaking (MMAMSB). This model is attractive theoretically because it can be derived 
from a concrete string compactification scenario due to Kachru, Kallosh, Linde and Trivedi 
(KKLT) which provides mechanisms for breaking SUSY, stabilizing unwanted moduli, and 
ensuring a positive cosmological constant [28]. In fact in Ref. [29] it is argued that the 
MMAMSB breaking pattern is a generic feature of a certain class of mechanisms that 
stabilize string moduli, of which the KKLT construction is only one example. 

A couple of papers have analyzed the soft SUSY breaking terms that arise from spe- 
cific realizations of this scenario, pointing out that the SUSY-breaking contributions from 
gravity (i.e. modulus) mediation are comparable to those coming from anomaly mediation 
[30, 31]. A parameter of the model, a, interpolates between pure anomaly mediation and 
pure modulus mediation. For vanishing a, the model suffers from the tachyonic slepton 
problem of pure AMSB, whereas for intermediate or large values of a, this problem is 
alleviated by the gravity-mediated contributions. 

The phenomenology of MMAMSB has been studied in several papers [32-39], where it 
has been shown that there are good prospects for discovering this type of SUSY breaking 
scenario at the LHC and/or at direct dark matter detection experiments. At least for 
certain choices of the model parameters, it should also be possible to measure them and 
identify the gaugino mass pattern distinctive to this model. 

In this paper, we aim to quantify the extent to which the MMAMSB scenario can be 
explored at the LHC by performing a global fit to prospective LHC measurements in the 
MMAMSB parameter space. This way, we can evaluate the accuracy and precision of the 
measurement of the model parameters. By attempting to fit different SUSY-breaking sce- 
narios, the CMSSM and mAMSB, to the MMAMSB model, we can also evaluate whether 
LHC measurements are sufficient to distinguish MMAMSB from other SUSY-breaking sce- 
narios. We show that the inclusion of rates as inputs to the fit is necessary to enable 
accurate parameter measurements with relatively low luminosity, and that rates are espe- 
cially crucial in distinguishing MMAMSB from the CMSSM and mAMSB. 

The layout of this paper is as follows. We begin by giving a brief introduction to the 
MMAMSB model and, in particular, the phenomenology of the gaugino sector. In Sec. 3, 
we discuss choosing a particular benchmark point in the model that satisfies all of the 
current experimental and observational constraints. In Sec. 4, we discuss how we fit the 
model with hypothetical LHC data and introduce the different observables that we use. 
Finally, in Sec. 5 we present the results of fitting the MMAMSB to observables derived 
from our benchmark point. We also try to fit other SUSY breaking models (the CMSSM 
and mAMSB) to our benchmark to see how quickly these models can be excluded and 
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which observables are effective at performing this task. We conclude with a discussion in 
Sec. 6. 

2 Mixed moduli-anomaly mediated supersymmetry breaking 

In this section we briefly introduce the MMAMSB scenario. We adopt the notation of 
Ref. [35]. As mentioned in the introduction, the MMAMSB can be derived from the 
KKLT string construction [28], or indeed from a class of string models in which moduli 
stabilization and SUSY breaking are accomplished dynamically [29]. In these models, the 
soft SUSY breaking terms receive contributions from both modulus and anomaly mediation. 
The parameter a interpolates between the pure modulus and pure anomaly extremes. 

The soft terms also depend on the sets of parameters rii and ia, which correspond to 
powers in the Kahler potential for matter fields and the gauge kinetic functions respectively. 
In the KKLT construction, £a can take on the value or 1 depending on whether the 
corresponding gauge field is localized on a D3 or D7 brane, whereas rii can be 0, 1, or 1/2 
for a matter field localized on a D7 brane, a D3 brane, or a brane intersection. 

Here we concentrate on the scenario where ia = ^ for all gauge fields and where rij = n 
for all matter fields (including the Higgs fields). Having a common value for rii among the 
sfermions is motivated by fiavor constraints. Having a different rii for the sfermions and 
the Higgses is possible, but for simplicity we assume a common value here. With these 
assumptions, the soft SUSY breaking parameters — the gaugino mass parameters, trilinear 
couplings, and sfermion mass parameters, respectively — at the GUT scale are given by 



Ma = j^[a + bagl], (2.1) 

Aijk = ^[3(^ - 1)« + in + 7i + 7fc)], (2.2) 

ml= (^)'[(l-n)a2 + 4ae,-7,]. (2.3) 

Here, ^3/2 is the gravitino mass, ga is the gauge coupling and ba = (^,1,-3) the 1- 



loop /? function coefficient for the gauge group o. The anomalous dimensions 7^ and their 
logarithmic derivatives 7^ as well as the mixed anomaly-modulus contributions are given 
in Appendix A of Ref. [33] . 

In Eq. (2.1) the term proportional to a is the universal gravity mediation contribu- 
tion, whereas the second term is the anomaly mediation contribution. The ratio of gauge 
couplings at the TeV scale should be a combination of the expected ratio from gravity 
mediation. 

Ml : M2 : M3 ~ 1 : 2 : 6 , (2.4) 
and the expected ratio from anomaly mediation. 

Ml : M2 : Afs ~ 3.3 : 1 : 9 . (2.5) 
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For the MMAMSB, it is 



,1 



Ml : M2 : M3 ~ (a + 3.3) : (2a + 1) : (6a - 9). 



(2.6) 



This distinct ratio of coupUngs is a hahmark sign of this model, and suggests that measuring 
the weak-scale gaugino masses could enable one to determine a. 

In pure anomaly mediation, there is a well-known problem of tachyonic slepton masses, 
which in the standard mAMSB scenario is solved by adding an ad hoc common term to 
all the sfermion masses. In this model, for moderate values of a the tachyonic slepton 
problem is solved by the gravity mediation contributions, and this term does not need to 
be included. 



In order to carry out a fit in the MMAMSB parameter space, we first need to select a 
point in this space to serve as our model. For a given point in parameter space, then, we 
must compute the spectrum and check that various constraints are satisfied and desirable 
features are present. 

We use the ISASUGRA spectrum generator that is packaged with ISAJET [40] to compute 
the RGE running of the parameters and to compute the low-energy spectrum and decays. 
We then use a number of codes to check various constraints: the IsaTools package pro- 
vided with ISAJET and the programs Superlso Relic [41] and micrOMEGAs [42-44]. Each 
of these codes computes the relic density and various precision low energy and flavor con- 
straints, and IsaTools and micrOMEGAs also compute dark matter direct detection cross 
sections. We also check that the model has not been ruled out by Higgs searches using 
HiggsBounds [45]. 

We consider scenarios with 1 ^ a 10, that have roughly similar SUSY breaking 
contributions from each sector and are close to the original KKLT construction value of 



In order for a fit using standard LHC observables, especially the kinematic endpoints, 
to be feasible, the decay chain 



must be present with a sufficiently large branching ratio. This means that the q ^ X2 ^i^d 
X2 ~^ decays must have sizeable branching ratios and sufficiently large mass splittings 
to give rise to hard enough jets and leptons to pass experimental cuts. In addition, the 
squark masses must be sufficiently low to give reasonable squark production cross sections 
at the LHC, but not so low that they are clearly ruled out by early LHC searches. 

^In this paper we follow the notation of Ref. [35] by Baer et al., in particular a = aBacr- As pointed out 
in that reference, however, the reader should be warned that elsewhere in the literature, e.g. Refs. [29, 32] 

by Choi et al., a is defined such that achoi = — ttt^t r — - — • The coupling ratio as a function of a will 

of course be a different function of achoi- 



3 MMAMSB model point 



a 



= 5. 
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PcirciiiiGtGr 


Vcllu6 


ParcinietGr 


Va 1 1 1 p 


a 


4.8 


Ml 


460.6 




21x10^ 


M2 


556.0 


tan/3 


10 


Ms 


976.7 


sign(/x) 


+1 






n 


0.5 







Table 1. MMAMSB benchmark point. All masses in GeV. 

For simplicity, we choose a common value of the modular weights rij for all of the 
sfermion and Higgs fields. In order to satisfy the branching ratio and mass splitting re- 
quirements just mentioned, we find that in this case n = 1/2 is the only valid choice. 
There is some freedom in choosing tan/3 and sign(/x), so we adopt the SPSla [46] values 
tan/3 = 10 and sign(/i) = +1. 

Performing a scan over the remaining parameters a and ?Ti3/2, we find that the strongest 
constraint on finding a viable model is the measurement by the WMAP experiment [47] 
of the dark matter relic density. The computed relic density takes on a value within the 
observed error bars in only narrow strips of the a — m^i2 plane^. One such region is an 
almost-vertical strip at a ~ 4.8, with m3/2 ranging from 15 TeV upwards, where a mostly 
bino though somewhat mixed x? undergoes efficient annihilation into h A, bb, W~ and 
Z H. From this region we select the benchmark point with = 21 TeV. The model 

parameters and the weak-scale gaugino mass parameters are summarized in Tab. 1. 

This model satisfies all the observational constraints implemented in the codes men- 
tioned above, though the spin-independent direct detection cross section is right on the 
edge of the XENONIOO [48] limit. The value of a is close to the value preferred by the 
KKLT string scenario (ckkklt = 5), and is far enough from the gravity- mediated and 
anomaly-mediated limits to give distinct phenomenology. In particular, the gaugino mass 
ratio (at the weak scale) in this model is Mi : M2 : M3 ~ 1 : 1.2 : 2.2, which is quite 
distinct from the ratio 1:2:6 expected from gravity mediation, and 3.3 : 1 : 9 expected 
from anomaly mediation. The model also has the branching ratios and mass differences 
mentioned above that are required for the standard LHC observables. 

The full spectrum of the model is given Tab. 2, and is plotted in Fig. 1 alongside the 
SPSla spectrum for comparison. Benchmark planes and lines for further LHC studies of 
MMAMSB models have been proposed in Ref. [49] 

4 Fit procedure 

In order to perform the fits we use the Fittino program [26]. Fittino first requires an 
input file, in which the user specifies the set of observables and provides their (perhaps 

^ For a given parameter point, the calculated relic density can vary significantly depending on which 
code is used to compute it and on which spectrum calculator is used. Because the relic density depends 
on the spectrum, and thus the model parameters, this difference can be undone by a slight change in the 
parameters and thus does not affect our results. 
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Table 2. MMAMSB benchmark model. All masses in GeV. 
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Figure 1. The spectrum of our MMAMSB model (on the left) compared to the spectrum of SPSla 
(on the right). Note the difference in scales on the vertical axis. 



hypothetical) measured values and uncertainties. Here the user also specifies the super- 
symmetric model and the high or low scale parameters to be fitted along with their starting 
values. Fittino then efficiently samples the parameter space, finds the best fit point and 
maps out confidence regions {i.e., contours of the likelihood function) using Markov chain 
Monte Carlo. 

It is important to emphasize that in the fits we perform here, we only use LHC ob- 
servables as inputs. While it is possible with Fittino to also include in the fits the various 
low-energy, flavor, and astrophysical observables that were used in the previous section to 
pick a benchmark model, we focus instead on the power of the LHC alone to constrain and 
measure the model parameters. 

In order to compute the likelihood function at each point in parameter space visited 
by the Markov chain, Fittino first calls an external code to calculate the mass spectrum, 
decay widths, and branching ratios. We use the ISASUGRA spectrum calculator, since it 
includes the MMAMSB scenario. 

The predicted values of the chosen observables are then computed internally by Fittino 
or by an external code, and compared with the measured values to compute the likelihood. 
Originally, Fittino included among its LHC observables a large set of kinematical edges, 
as well as some branching ratio observables. Crucial to our analysis will be the inclusion 



-7- 



of event rates, which, as mentioned earher, were implemented as observables in Fittino 
in Ref. [25]. 

In this implementation, Fittino first uses the decay table computed by the spectrum 
calculator to determine the branching ratios of the squark- and gluino- initiated decay chains 
that could contribute to the given channel. Each branching ratio is then multiplied by the 
relevant squark and/or gluino production cross section. This cross section is interpolated 
from a grid over rriq and rrig of NLO cross sections that have been previously computed 
using the code Prospino [27, 50, 51]. It remains to compute the acceptance, i.e. the fraction 
of events which pass the analysis cuts for the given channel. As is described in detail 
in Ref. [25], in this implementation the acceptance is computed by combining analytical 
formulae for the momentum distributions of the decay products of at-rest squarks and 
gluinos with numerical estimates (which are also interpolated from pre-computed grids) of 
the effect of boosting the system into the lab frame. 

It was shown in Ref. [25] that this fast technique for estimating the acceptance agrees to 
within 5% with acceptances determined using Monte Carlo event generation. To verify that 
this agreement holds within the parameter space of the MMAMSB scenario, we computed 
acceptances across a grid in MMAMSB parameter space using Herwig++ [52, 53] to generate 
events and Rivet [54], including the anti-fct jet finder [55, 56], to apply the cuts. We then 
performed some of our MMAMSB fits using acceptances interpolated from these grids 
alongside the fits using the rates code implemented in Fittino, and found very good 
agreement. 

An extremely important ingredient in the calculation of the likelihood is the estimate 
of the uncertainty on each observable. For the well-studied SPSla model, which is the basis 
of the Fittino fits performed in Refs. [19, 25], the uncertainties for most observables can 
be obtained from the thorough study carried out in Ref. [17]. For the fit we perform here, 
we must extrapolate these uncertainties from SPSla to our MMAMSB model. To do this, 
we must take into account the relative number of events from SPSla and our model in each 
signal channel. In our MMAMSB model, the squarks and gluino have higher masses than 
they do in SPSla, so the production cross sections are lower, leading to fewer signal events 
and thus larger uncertainties on most observables. In some cases, however, an increase in 
branching ratios between models can more than make up for the decrease in cross section. 
For example, the branching ratio of the decay xl ~^ ten times larger in our model 

than in SPSla, leading to a much more pronounced di-lepton edge due to xt decays. 

To verify that the scaling of uncertainties we performed is reasonable, we looked at one 
of the observables, the mu endpoint (for a more detailed description of each observable, 
see the next two subsections) in detail. We generated events for SPSla and our MMAMSB 
model, including energy smearing to account for detector effects. We then fit a smeared 
step function to the data in order to extract the endpoint and determine its uncertainty. 
The uncertainties we obtained in this way roughly match the SPSla uncertainties given 
in Ref. [17] and the uncertainties we extrapolated from this for our MMAMSB model, 
verifying the reasonableness of our extrapolations. For illustration, in Fig. 2 we show the 
mu distribution, with statistical error bars, for each model for a particular choice of center- 
of-mass energy and integrated luminosity, with the fitted function superimposed over the 
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SPSIA, 14 TeV, 10 fb~' MMAMSB, 14 TeV, 10 flj'' 




m„ (GeV) m„ (GeV) 



Figure 2. The distribution for SPSla (left) and the MMAMSB model (right) for 10 fb^i of 
integrated luminosity at the LHC with \/s — 14 TeV. The error bars show the statistical error only. 
Superimposed is the fit function used to extract the endpoint value. 

data. 

4.1 Kinematic edges 

Throughout this paper, we compare the impact of various sets of observables on the fits. 
The first two sets below do not include rates, and are mostly edges of kinematical distribu- 
tions. The basic set we call Group I and includes standard kinematical edges built from the 
"golden" decay chain described earlier. Group II contains additional quantities built from 
this decay chain as well as some observables sensitive to the properties of third generation 
sparticles. The full list of Group I and II observables is as follows: 

• List of observables in Group I (these are all defined in Ref. [16]). 

— m^^^, the dilepton invariant mass edge. 

— m^^f^, the jet dilepton invariant mass edge. 

— m^gf", the jet-lepton low invariant mass edge. 

— iTi^'^^, the jet-lepton high invariant mass edge. 

• List of observables in Group II. 

— ?7T-*^£, the jet-dilepton threshold invariant mass edge [16]. 

— m?^, the squark stransverse mass [57, 58]. 

— m^^^, the di-tau invariant mass edge [17, 19]. 

— mjj, the weighted top-bottom invariant mass edge [19]. 

— Anig^o, the mass difference between the gluino and the LSP [17, 19]. 

— '^(^o^£^; the dilepton invariant mass edge from the decay of a Xa [17, 19]. 

— r?^, the ratio of selectron (smuon) to stau mediated X2 decays [19]. 
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It 


0.457 
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0.0036 







Table 3. LHC observables for the MMAMSB benchmark point, Tab. 1. The masses and branching 
ratios have been calculated with ISASUGRA [40]. The uncertainty estimates on the observables are 
based on [17, 19] and have been rescaled as described in the main text. All dimcnsionful quantities 
are given in GeV. 

In Tab. 3 we provide, for each observable, its nominal value in our MMAMSB model 
and its uncertainty, determined by scaling from SPSla as described above, for each luminos- 
ity and center-of-mass energy scenario. The last two columns give additional contributions 
to the uncertainty from the lepton energy scale and jet energy scale [19]. 

4.2 Rates 

In addition to the Group I and II observables listed in the previous subsection, we also 
consider the two rate observables discussed previously. Here, we give the description of 
these observables, including the full list of cuts. 

• List of rate observables. (These observables are also defined in Ref. [25].) 

— Rjj^^, the event rate for at least two hard jets and missing transverse energy. 

* PTjct > 50 GeV. 

* |7?jet| < 2.5. 

* i?jet = 0.4 (anti-A:^ jet algorithm [55, 56]). 

* > 100 GeV. 

— R^ijjf^, the event rate for at least two hard jets and missing transverse energy 
with a pair of opposite sign, same flavour leptons (1st or 2nd generation). The 
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Observable 


7TeV 

Value (fb) Uncertainty 


14 TeV 
Value (fb) Uncertainty 




113 23 
11.8 3.5 


2780 556 
245 49 



Table 4. LHC event rates for the MMAMSB benchmark point, Tab. 1. The event rate includes 
the NLO squark and gluino production cross section [50, 51], the branching ratios of the decays 
and the expected particle acceptances. The acceptances were tested with full parton shower and 
hadronisation using Herwig++ [52, 53], Rivet [54] and the anti-fct jet finder [55, 56]. 

background from the leptonic decays of r leptons, Xi^2 ^^"^ removed by 

subtracting events with opposite sign, different flavour lepton pairs. 

* ^jjtr signal. 

* PT,e > lOGeV. 

* [%[ < 2.5. 

* Lepton-jet isolation of AR = 0.2 (jet activity within AR = 0.2 < 10 GeV). 

In Tab. 4, we provide the values for these observables, i.e. the event rates for the two 
different types of events as predicted by a Monte Carlo simulation. We provide the rates 
and their uncertainties for both 7 and 14 TeV. 

5 Fit results 

In this section we present the main results of our paper and show that including event 
rates in fits can significantly increase their effectiveness to measure and exclude different 
SUSY breaking scenarios. We begin by fitting to our own MMAMSB benchmark scenario 
and show that even with early data, the model can be constrained and some parameters 
measured accurately. With larger data sets, we should be able to measure the whole of the 
parameter space with high accuracy (< 5%). 

We then try to fit other SUSY breaking scenarios to our benchmark point from 
MMAMSB. We start with the CMSSM and show that if we only include the Group I 
mass edges, a good fit can be achieved unless we have very large data sets (100 fb~^ at 
14 TeV). However, if we include event rates in the fit, we see that even with a small amount 
of data (10 fb"^ at 7 TeV), we can exclude the CMSSM. We complete the same task for the 
mAMSB and see very similar results. With Group I observables and early data we cannot 
conclusively exclude the scenario. However, as soon as we add event rates, the model can 
no longer fit the data. 

5.1 MMAMSB 

In this section we discuss the results of fitting the MMAMSB model to observables derived 
from our benchmark scenario. Fits are performed for 10 fb'^ at 7 TeV and 1 fb-\ 10 fb-\ 
and 100 fb"^ at 14 TeV. In addition, we also show the effect of adding rates (see Tab. 4) 
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Figure 3. 1-a (68.3%) (black), 2-cr (95.4%) (grey), and 3-a (99.7%) (light grey) two-dimensional 
confidence regions for 7 TeV, 10 fb~^, using Group I observables, without (left column) and with 
(right column) rates. The green 'X' represents the input point and the red circle the best-fit point. 



to the fits, and including the more complicated set of Group II observables (see Tab. 3). 
We demonstrate that rates significantly improve the accuracy of the fits, especiahy when 
using early LHC data with limited statistics. 

Let us first consider the fits that can be performed with 10 fb~^ at 7 TeV and 1 fb~^ 
at 14 TeV. With this amount of data and our benchmark scenario, we do not expect to 
have sufficient statistics to be able to measure any of the observables given in Group II. 
Therefore, we perform the fits using the mass edges in Group I, both with and without the 
rates. In the left-hand column of Figs. 3 and 4, we see that when only using the Group 
I edges, the model is essentially unconstrained across the whole parameter space and an 
effective fit cannot be performed. 

As soon as we add rates to the fit, however, the situation improves remarkably. At 
our benchmark point we can now constrain 171^/2 to ~ 15% and a to ~ 10% as can be seen 
in the right-hand column of Figs. 3 and 4, or in the zoomed-in version of these regions 
shown in Fig. 5. The fact that the rates improve the constraints on ^71,3/2 and a is to 
be expected if we examine the form of the soft breaking terms in the MMAMSB model, 
Eqs. (2.1)-(2.3). We see that the soft breaking masses are all proportional to the product 
Therefore, the rate observables, which are sensitive to the overall mass scale of 
any new states, constrain this combination effectively. The modular weight, n, is less well 
constrained — at 7 TeV it is only determined within ~ 20%, as can be seen in Fig. 5 — due 
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Figure 4. Like Figure 3 except for 14 TeV, 1 fb ^. 



to the fact that it has a sub-leading dependence in the soft breaking terms. This is stih 
an improvement, however, over the fit without rates, where n could not be constrained 
at all. Finally, tan/3 has the worst constraints of all the parameters in the fit; even with 
rates included, most values are allowed. As we show later, a measurement of tan /3 requires 
observables that are sensitive to third-generation particles and none of the observables in 
Group I measure these. 

If we move to the fits with 10 fb^^ at 14 TeV, we see in Fig. 6 that even an accurate 
determination of all the Group I mass edges can still not constrain the MMAMSB model 
in any parameter. Once we add the information from rates, we again see a good fit for 
m3/2, ct and n, and even tan/3 is constrained (tan/3 



10l'^3'5). "^^^ reason that tan/3 can 



be constrained with the rate observables is that the combination of R 



and Reejjf^ acts 



as an observable for the branching ratio of the decay X2 ~^ This branching ratio 

is sensitive to tan/3 since it competes with the branching ratio of the decay X2 ~^ t^t^, 
which increases substantially as tan f3 becomes large. 

Adding the Group II set of observables with a reasonable statistical error becomes 
possible with 10 fb~^ at 14 TeV. We can now compare the fit using the Group I and II 
observables alone with that of using just the Group I observables along with rates. We 
see in Fig. 7 that in general the fit is greatly improved, especially for tan/3 and n. The 
constraining power for tan/3 mainly comes from the ratio r?3 and the mass edge m,-^. 
These observables directly measure the contribution of third generation superpartners, 
which is most sensitive to tan /3. The measurement of the mass scale is not improved 
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Figure 5. A zoomed-in view of the confidence regions for 7 TeV, 10 fb ^ and 14 TeV, 1 fb 
These are the same fits as in Fig. 3 and Fig. 4. 



significantly, however, as this was already well constrained by the rate information. In the 
Group II observables, the squark stransverse mass {m'^'^) and the jet-dilepton threshold 
{iTi^^[) now perform this task. 

Despite the extra observables that constrain the mass scale in Group II, if we add rates 
to the fit, the measurement of both improves substantially, as is clear from 

Fig. 7. As stated before, this is due to the rates being particularly sensitive to the overall 
mass scale. The measurement of 171^/2 is improved by almost a factor of 2. 

An interesting feature of the fits that now becomes much more clear is the double 
minimum seen in the tan/3-n plane. The double minimum is due to the fact that the 
functional forms for m^f^, m^gi" and rn^^^^ depend on the relative mass ratios of the 
particles in the cascade decay [15, 16]. For example, in the correct minimum we have a 
mass ordering 2m? > m-o+m^o > 2m^om^o, which leads to the end-points given in Tab. 3. 
However, in the second minimum, the mass ordering switches to mln +W'?o > 2m^om^o > 
2m|, where the end-points have a different functional form. Thus, a different set of masses 
can lead to the same end-points being measured at the LHC. The problem of end-point 
mimics was previously been discussed in Refs. [59, 60] and possible solutions proposed 
were to measure the Xi mass accurately at a linear collider [61, 62] or to measure the whole 
invariant mass distribution as the shape would be different for the mimic distribution. 
Unfortunately, the rates are unable to resolve the end-point mimics as they are only weakly 
sensitive to the masses of the electroweak particles further down the decay chain. The rates 
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Figure 6. Confidence regions for 14 TeV, 10 fb ^, only including Group I observables, with (left 
column) and without (right column) rates. 

are most sensitive to the first generation squarks and gluino masses and in both minima, 
these are very similar. 

With 100 fb-i of data, the Group I mass edges are all measured with high accuracy. 
Consequently, these mass edges alone can constrain all of the parameters apart from tan /3 
to some extent, as can be seen in Fig. 8. Adding the rates to the fit, though, improves 
the measurements significantly. We see a large improvement in the determination of 
(the width of the la confidence region drops from 10 TeV to 4 TeV) and a, while tan/3 is 
becoming constrained. 

If we include Group II observables to the above fits, as shown in Fig. 9, we now 
achieve a precision determination of better than 5% on all parameters. In fact, a, which 
parameterises the ratio of gauge couplings [see Eq. (2.6)] is measured to better than 1%. In 
addition, the double minimum displayed in other fits is no longer present and we now have 
no mass ordering ambiguity. The observable that breaks the ambiguity is m°^^ = 32.0 GeV 
(measured 58.9 GeV) which is 5.3 a away. This is tensioned against m^^^ which is the most 
accurately measured observable. If we want to increase m^^^ for this point, we would also 
increase m^^^, leading to an even worse fit. 

In the high accuracy fits using 100 fb~^ (Figs. 8 and 9), two interesting correlations 
become apparent. First, in the m3/2~a plane, there is a negative correlation between the 
two parameters. This correlation can be easily understood when we inspect the form of 
the soft breaking terms in the model, Eqs. (2.1)-(2.3). All of the soft masses are oc 
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Figure 7. Like Fig. 6, but now also including Group II observables. 
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Figure 8. Like Fig. 6, but using 100 fb ^. 
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Figure 9. Like Fig. 7, but using 100 fb ^. 



Therefore to keep the particle masses at a set value, as ^3/2 increases, a must decrease 
and vice versa. 

The other clear correlation is a positive one between the parameters tan (3 and n. This 
correlation is due to the trilinear coupling Aijf. which enters in the f mixing matrix. The 
mixing term in the f sector is 



m- 



mT-{Ar — /xtan f3) 



(5.1) 



and from Eq. (2.2) we can see that Ar oc n. Therefore, in order to keep the mixing in the f 
sector constant, if n {Aij^) increases tan/3 must also increase to compensate. The opposite 
would be true of the i sector but the observables r?^ and m™^^^ are far more accurately 
measured than the observable mf^. Thus the fit is dominated by the constraint in the f 
sector. 

The results of our MMAMSB fit are collected in Tab. 5. 
5.2 CMSSM 

In the previous section we have shown that we can accurately reconstruct our MMAMSB 
benchmark point at the LHC, even with relatively small amounts of data (10 fb^^ at 
14 TeV). Here we would like to determine, however, if it is also possible to rule out other 
SUSY breaking scenarios. Possibly the most widely studied scenario is the CMSSM and 
thus it is natural to ask what happens if we attempt to fit this scenario to our benchmark 
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Table 5. Fits to MMAMSB parameters for our chosen benchmark point. Fits are done with various 
sets of observable groups (I and II) and errors (Tab. 3). Fits are also done with and without the 
rates observables (Tab. 4). 



point. In particular we would like to use the of the best-fit point to determine if the 
CMSSM can be ruled out. 

We begin by performing the fit for 10 fb^^ at 7 TeV using only the Group I observables 
and without the information from rates. ^ The results of this fit are shown in Tab. 6. We 
see that the fit has a viable minimum with a of just 0.12. In addition, the fit begins 

^For these fits the number of measurements is equal to the number of free parameters. Therefore, a 
perfect fit to the data (x^ ~ 0) should in principle always be possible, but this is not the case, as is shown 
by the fits at 14 TeV with 10 fb"^ and 100 fb"^ where ^ = 2.1 and = H-S respectively There are 
a couple of reasons for this. One is that the Group I observables only measure the wino, bino and first 
generation sfermion masses. Thus, they only have weak sensitivity to tan /3 and ylo and the fit is effectively 
to only two free parameters. A second reason is that many other constraints are implicitly present in the 
fit but not formally included in the number of measurements. For example, we require a neutral LSP so 
that we have a dark matter candidate, consistent electroweak symmetry breaking, and no tachyonic degrees 
of freedom. All of these constrain the areas of parameter space that we are able to search to satisfy the 
mass-edge conditions. 
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Table 6. Best fit points for the CMSSM and tlie minimum for that point and associated set of 
observables. We only include the 1-a environment when the best fit point is not excluded at the 
99.9% confidence level. We see that rates are extremely effective at ruling out the CMSSM. 



to constrain the parameter mi/2 (160 < mi/2 < 220) while the other parameters can vary 
quite freely, as is illustrated in Fig. 10. As we move to 1 fb~i at 14 TeV but keep the 
same observable set, constraints begin to appear on ttt-q (0 < rriQ < 100) and we still have 
a minimum where all the observables are wen fitted (x^ = 0.23). 

If we add the rate information to the fit, however, the result changes spectacularly 
and the best fit has = 216 which means that the point is completely excluded by the 
data. If we examine the individual measurements at this best fit point more closely, we 
see that the exclusion is dominated by two measurements. The jet-lepton high invariant 
mass edge, whose true value is tti^^^^ = 652 GeV, has a measured value of 311.9 GeV at 
the best-fit point, a full 13 a away. The event rate at the best fit point, iJ^-^^miss = 231 fb 
(measured = 113 fb), is 5.3 a away. Essentially, the CMSSM model cannot replicate the 
more compressed mass spectrum with the relatively large masses for the gluino and squarks. 
The edge observables (especially the difference between the squark and slepton mass) are 
trying to pull the mass scale down, while the rate observables try to pull the mass scale 
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Group I 




Figure 10. Plots showing the fit of the CMSSM to our MMAMSB benchmark point using only 
Group I observables and no rate information. The fits in the upper row are for 7 TeV, 10 fb~^ while 
the bottom row shows the fits for 14 TcV, 1 fb^^. 

up. The tension between these two observables leads to the exclusion. 

We can also examine the best-fit point found using the Group I mass edges alone to see 
what rates are predicted for these points. We find that for the best-fit point with 10 fb~^ 
at 7 TeV, Rjj^mms = 14809 (measured = 113) and is 652 a away! Incidentally, this best fit 
point is already easily ruled out by current searches [5, 7, 12, 13]. 

Moving to the fits performed at 14 TeV with 10 fb~^, we see that using the Group I 
observables alone, = 2.1, and thus the fit still has a viable minimum. If we add the 
Group II observables, however, the best-fit point now has a x^/d.o.f.= 122/7 and is thus 
excluded. The exclusion is dominated by the stransverse mass observable, which at the 
best-fit point is m?^ = 444 GeV, 8.9 a away from the measured value m'^'^ = 699 GeV. 
The stranverse mass is effective as its rough dependence is mj^ oc y/m? — 2m|o and thus 
with heavy squarks, it is essentially measuring the squark mass. This leads to a similar 
tension as in the fits with rates: the stransverse mass tries to pull the mass scale of the fit 
up while the mass edges try to bring the scale down. 

Even though the exclusion due to the stransverse mass seems convincing, we should 
compare this to the exclusion that would be achieved by using the rate observables. At the 
best-fit point, the inclusive rate has the predicted value -Rj^gmiss = 125 pb, which is a huge 
218 a away from the value Rjj^miss = 2.78 pb that would be measured. Thus we can see 
that the rates are far more effective observables for determining the overall mass scale of a 
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Table 7. Best fit points for the mAMSB and the minimum for that point and associated set of 
observables. We only include the 1-a environment when the best fit point is not excluded at the 
99.9% confidence level. We see that rates are extremely effective at ruling out the mAMSB. 

model than the stransverse mass is. This can also be seen by comparing the fit done with 
10 fb~^ at 7 TeV with only the Group I observables and including rates with the fit done 
with 10 fb~^ at 14 TeV with Group I and II observables but no rate information. From 
Tab. 6, we see that even though the fit with rates is done with fewer observables in total 
and far less data, with a total = 216 it more convincingly rules out the CMSSM model 
than the Group II fit (total = 122) that would take at least a few years longer. 

5.3 mAMSB 

In the previous section we have shown that the MMAMSB can be convincingly distin- 
guished from the CMSSM even with early LHC data (10 fb"^ at 7 TeV). We now do a 
similar analysis witha different SUSY breaking scenario, mAMSB, to see if the LHC can 
perform the same task and also separate this model. The mAMSB scenario nearly corre- 
sponds to pure anomaly mediation, i.e. the q — )• limit of the MMAMSB. As mentioned 
earlier, however, pure anomaly mediation leads to tachyonic slepton masses. In minimal 
AMSB (mAMSB), this problem is solved with the addition of a constant contribution mo 
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Figure 11. Plots showing the fit of mAMSB to our MMAMSB benchmark point using only Group 
I observables and no rate information. The fits in the upper row are for 7 TeV, 10 fb~^ while the 
bottom row shows the fits for 14 TeV, 1 fb^^. 



to the scalar masses at the GUT scale. Since a = means that the n-dependence in the 
model also disappears, mAMSB has only three parameters, mo, ^,3/2, and tan/3. 

We begin by performing the fit with only the Group I observables at 10 fb~^ at 7 TeV. 
The results of this fit are shown in the upper row of Fig. 11. The fit has a viable minimum, 
but a x^/d-o.f. = 3.8/1. Even with just these four mass edge observables and early data, 
the model already suffers from a little tension. In addition, both mo (100 < mo < 150 GeV) 
and m3/2 (13.4 < m3/2 < 16.4 TeV) are already relatively constrained. As we move to 
1 fb~^ at 14 TeV, the tension in the fit increases and the minimum now has a x^/d-o.f. 
= 7.6/1. This corresponds to a p-value of less than 0.01 and shows the model is already 
highly disfavoured. The reason why the Group I observables alone can exclude the model is 
that the mAMSB contains a larger mass splitting between the gaugino masses than either 
MMAMSB or the CMSSM (see Sec. 3). 

If we include the rates in the fit, however, the exclusion of the model becomes far 
clearer, as can be seen in the lower row of Fig. 11. Including the event rates with 10 fb~^ 
at 7 TeV leads to a best fit point with a minimum x^/d.o.f. = 238/3 and the model is 
completely ruled out. The reason is the same as for the CMSSM: to fit the mass edges 
requires the mass scale in mAMSB to be low. For example, the masses of all the coloured 
particles in the best fit point of the Group I fit with 10 fb~^ at 7 TeV are under 400 GeV 
and the event rates are far higher than they are for the model we are fitting to. This 
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analysis is confirmed when we look at the individual observables. At the best fit point the 
jet-lepton high invariant mass edge has the value fn^^^ = 653 GeV, which is 13 a away 
from its value (312 GeV) at the MMAMSB benchmark point, and is trying pull the mass 
scale down. On the other hand, the event rate at the best fit point, Rjj^miss = 240 fb, is 
5.7 a away from the MMAMSB value of 113 fb and is trying to pull the mass scale up. 
The tension between the two observables leads to the clear exclusion. 

As for the fit to the CMSSM, we can also examine the effect of the Group II observables 
on the exclusion and in particular the squark stransverse mass, MJ'^, which is sensitive to 
the mass scale. We find the fit with Group I and II observables but without event rates 
with 10 fb^"*^ at 14 TeV has x^/d.o.f. = 330/8 and thus excludes the model. As expected, 
the exclusion is dominated by MJ'^ which is 12 o" away from the true value at the best fit 
point. Once again, we have a tension between the edge observables, which are trying to 
pull the mass scale down, and MJ^, which is trying to pull the mass scale up. 

In the CMSSM fit, we also compared the effectiveness of the squark stransverse mass 
observable with that of the event rates. Here we do the same by taking the above best fit 
point with just Group I and II observables and examining what the event rate would be 
at this point. We find that the event rate Rj^^miss = 539 pb for the best fit point which 
is 953 a away from the value (2.78 pb) in our MMAMSB benchmark scenario. Again, this 
confirms that rates are far more sensitive to the mass scale than rn^"^ and provide a much 
more constraining test of the SUSY breaking scenario. 

6 Discussion 

In the MMAMSB both modulus and anomaly mediated terms can contribute at roughly 
equal levels to SUSY breaking. This model can have very similar phenomenology to the 
widely studied CMSSM and mAMSB models. A unique feature of the MMAMSB, however, 
is that depending on the value of the phenomenological parameter a, the ratio of gaugino 
masses can be very different from other SUSY breaking scenarios. 

In this paper we determined how well we can expect the MMAMSB model to be mea- 
sured at the LHC and if it is possible to distinguish it from other SUSY breaking scenarios. 
We fit the model to various hypothetical LHC measurements to see how precisely the in- 
put parameters can be determined. Among the inputs to the fits, we included the widely 
studied mass edges of various SUSY cascade decay chains, as well as other kinematical 
observables and measurements of ratios of branching ratios that are expected to be pos- 
sible at the LHC. Special to our procedure, we also included event rate observables, to 
try and further constrain the model parameters. Event rates are extremely sensitive to 
the mass scale of the model but are conventionally difficult to include in fits due to the 
computationally intensive task of running a Monte Carlo event generator for each point 
in the fit. This is solved by using a combination of cross-section grids and analytically 
calculated acceptances that is vastly quicker [25] . 

Using a particular MMAMSB benchmark scenario we showed that the parameters of 
the model will be difficult to reconstruct with early data if only mass edge measurements 
are used. When including the event-rate information, however, even early LHC data is 
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sufficient to constrain the model, especially the parameters that determine the overall 
mass scale. With more data and more intricate observables, we showed that all the model 
parameters can be reconstructed with high accuracy (< 5%). 

In addition to showing that the MMAMSB model can be reconstructed, we also demon- 
strated that it can be distinguished from the CMSSM and the mAMSB. Using only the 
basic 'Group V mass edges (see Sec. 4), we saw that even with 100 fb~^, it is difficult to 
conclusively exclude the CMSSM from being able to fit our MMAMSB benchmark point. 
With the addition of event-rate information, however, the CMSSM interpretation is com- 
prehensively excluded with just 10 fb~^ at 7 TeV. We reached a similar conclusion when 
trying to fit the mAMSB to our benchmark point. Again, when we included rates, the 
model was easily excluded with early LHC data. 

We have seen that rates play a crucial role in measuring and distinguishing different 
SUSY models. The models considered in this paper have few free parameters. Therefore 
an obvious direction for future work is to investigate more general models and see if rate 
information can allow the precise determination of some or all of their parameters. 
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